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\ Abstract 

> , 

Q For d > 1 and a G (0, 2), consider the family of pseudo differential operators {A + 6A Q / 2 ; b 6 

^ | | [0,1]} on M. d that evolves continuously from A to A + A Q / 2 . In this paper, we establish a 

uniform boundary Harnack principle (BHP) with explicit boundary decay rate for nonnegative 
functions which are harmonic with respect to A + 6A"/ 2 (or equivalently, the sum of a Brownian 
motion and an independent symmetric a-stable process with constant multiple in C 1,1 open 
sets. Here a "uniform" BHP means that the comparing constant in the BHP is independent 
p | of b G [0, 1]. Along the way, a uniform Carlcson type estimate is established for nonnegative 

functions which are harmonic with respect to A + 6A Q / 2 in Lipschitz open sets. Our method 
employs a combination of probabilistic and analytic techniques. 
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1 Introduction 

> ■ 

^ 1 Discontinuous Markov processes have been receiving intensive study recently due to their impor- 

tance both in theory and in applications. Many physical and economic systems could be and 
in fact have been successfully modeled by discontinuous Markov processes (or jump diffusions as 
some authors call them); see for example, |28t [33| [35] and the references therein. The infinites- 
imal generator of a discontinuous Markov process in M. d is no longer a differential operator but 
rather a non-local (or integro-differential) operator. For instance, the infinitesimal generator of a 
rotationally symmetric a-stable process in M. d with a G (0, 2) is a fractional Laplacian operator 
cA a / 2 := -c(-A) Q / 2 . 

Discontinuous Markov processes include the very important Levy processes as special cases 
and they are of intrinsic importance in probability theory. Integro-differential operators are very 
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important in the theory of partial differential equations. Most of the recent study concentrates 
on discontinuous Markov processes, like the rotationally symmetric a-stable processes, that do not 
have a diffusion component. For a summary of some of these recent results from the probability 
literature, one can see \15\ [TO] and the references therein. We refer the readers to [12\ \T3[ Q3] for a 
sample of recent progresses in the PDE literature. 

However, in many situations, like in finance and control theory, one needs Markov processes 
that have both a diffusion component and a jump component, see for instance, [27 ^1341135] , The fact 
that such a process X has both diffusion and jump components is the source of many difficulties in 
investigating the potential theory of the process X. The main difficulty in studying X stems from 
the fact that it runs on two different scales: on the small scale the diffusion part dominates, while 
on the large scale the jumps take over. Another difficulty is encountered when looking at the exit 
of X from an open set: for diffusions, the exit is through the boundary, while for the pure jump 
processes, typically the exit happens by jumping out from the open set. For the process X, both 
cases will occur which makes the process X much more difficult to study. 

Despite these difficulties, in the last few years significant progress has been made in understand- 
ing the potential theory of such processes. Green function estimates (for the whole space) and the 
Harnack inequality for a class of processes with both continuous and jump components were estab- 
lished in (36] and [37]. The parabolic Harnack inequality and heat kernel estimates were studied 
in [38] for Levy processes on M rf that are the independent sum of Brownian motion and symmetric 
stable process, and in [20] for much more general symmetric diffusions with jumps. Moreover, a 
priori Holder estimate is established in [20] for bounded parabolic functions. For earlier results on 
second order integro-differential operators, one can see [23] and the references therein. 

The boundary Harnack principle (BHP) is a result about the ratio of positive harmonic func- 
tions. We say that the BHP holds for an open set D C M. d if there exist positive constants i?o and 
C depending on D with the property that for any Q E dD, r £ (0, Ro], and any positive harmonic 
functions u and v in D n B(Q, r) that vanish continuously on dD n B(Q, r), we have 

^l<C^rl for allx,y £ DnB(Q,r/2). (f.f) 
u{y) v{y) 

The BHP for Brownian motion (or, equivalently, for the Laplacian) is a fundamental result 
in analysis and PDE. It was independently established for Lipschitz domains in the late 1970's by 
Ancona, Dahlberg and Wu ( [H [221 132] ) • Later, Bass and Burdzy developed a probabilistic method in 
[5] to prove the boundary Harnack principle and extended the boundary Harnack principle to more 
general domains (see also [4]). When D is a bounded C 1 ' 1 domain, (jl.lj) can be strengthened to the 
following version that gives the explicit boundary decay rate of non-negative harmonic functions 
that vanish on the boundary: 

44<Ct44 ioraHx,y€DnB(Q,r/2), (1.2) 

u{y) °D{y) 

where 5d{x) is the Euclidean distance between x and D c . The BHP plays a vital role in the study 
of potential theory of Brownian motion and Dirichlet Laplacian in domains. For example, BHP 
can be used to show that Martin boundary can be identified with the Euclidean boundary for a 
large class of domains and to study the non-tangential limit of non-negative harmonic functions 
near the boundary (see [2] for an analytic approach and [3] for a probabilistic approach). In fact, 
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BHP has also be established for a large class of diffusions (or, equivalently, for second order elliptic 
equations), see (TT| [23"]. 

The study of BHP for discontinuous Markov processes and integro-differential operators is quite 
recent. It was first established for bounded Lipschitz domains in [7] and then extended to more 
general open sets in [40] . Subsequently Bogdan-Stos-Sztonyk [9] and Sztonyk [4T] extended the 
boundary Harnack principle to symmetric (but not necessarily rotationally invariant) stable pro- 
cesses. Recently, the BHP has been extended in [31] to a large class of pure jump Levy processes 
that can be obtained from Brownian motion through subordination. Very recently, the boundary 
Harnack principle for some one-dimensional Levy processes with both continuous and jump com- 
ponents was studied in [32]. However BHP for processes on R d in dimension two and higher that 
have both diffusion and jump components have been completely open until now. Note that the fact 
that a pure jump process may (and typically does) exit an open set by jumping out of it stipulates 
that, in the boundary Harnack principle for such processes, the nonnegative harmonic functions 
vanish continuously on D c n B(Q, r). 

The principal goal of this paper is to establish the boundary Harnack principle for nonnegative 
functions which are harmonic with respect to the independent sum of a Brownian motion and a 
symmetric stable process in C 1,1 open sets in M. d for every d > 1. The process X studied in this 
paper, although quite specific, serves as a test case for more general processes with both continuous 
and jump parts. The study of this test case will hopefully shed new light on the understanding of 
the boundary behavior of nonnegative harmonic functions of general Markov processes. 

Intuitively, the independent sum X of a Brownian motion and a symmetric stable process can 
be thought roughly as some sort of "perturbation" of Brownian motion. Thus some people might 
expect the BHP for X could be established through some general perturbation technique. However, 
this kind of approach may not always work. In [291 [30] , the potential theory of truncated symmetric 
stable processes including BHP was studied. One of the main results in [29] is that the BHP is valid 
for the positive harmonic functions of this process in bounded convex domains. A very interesting 
fact is, even though truncated symmetric stable processes can be considered as a perturbation of 
rotationally symmetric stable processes (see |25[ 130]). unlike symmetric stable processes, the BHP 
for truncated symmetric stable processes fails in non-convex domains (see the last section of |29] 
for a counterexample). This indicates that general perturbation method may not be suitable for 
establishing the BHP. 

Let us now describe the main result of this paper more precisely and at the same time fix the 
notations. A (rotationally) symmetric a-stable process Y = (Yt,t > 0,¥ x ,x £ M. d ) in ]R rf is a Levy 
process such that 



The infinitesimal generator of a symmetric a-stable process Y in M. d is the fractional Laplacian 



where A(d, a) := a2 a " 1 7r~ ci / 2 r(^)r(l - § Here T is the Gamma function defined by T(A) := 
/ °° t x ~ l e^dt for every A > 0. 



E x e iHYt ~ Yo) = e~ ma for every x G R d and £ <E R d . 



A"/ 2 , which is a prototype of nonlocal operators. The fractional Laplacian can be written in the 
form 
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Suppose X° is a Brownian motion in R d with generator A = Y2i=i ~jh?-i anc ^ ^ is a symmetric 

i 

a-stable process in R d . Both X° and Y satisfy a self-similarity, which will be used several times in 
this paper. That is, for every A > 0, {A^ 1 / 2 ^ - Xg), t > 0} and {^^{Yxt - Y ), t > 0} have the 
same distributions as that of {X t ° — X§,t > 0} and {Y t — Yo,t> 0}, respectively. Assume that X° 
and Y are independent. For any a > 0, we define X a by Xf := X® + aYt. We will call the process 
X a the independent sum of the Brownian motion X° and the symmetric a-stable process Y with 
weight a > 0. The infinitesimal generator of X a is A + a a A a / 2 . For every open subset D C R rf , we 
denote by X a,D the subprocess of X a killed upon leaving D. The infinitesimal generator of X a,D 
is (A + a a A a ^ 2 )\u- It is known (see [38]) that X a > D has a continuous transition density p%(t, x, y) 
with respect to the Lebesgue measure. We will use p a (t,x,y) to denote the transition density of 
X a (or equivalently, the heat kernel of A + a a A a l 2 ). The quadratic form (S,J-) associated with 
the generator A + a a A a / 2 of X a is given by 

T = W 1 ' 2 (R d ) := \ u G L 2 (R d ;dx) : ^- G L 2 (R d ;dx) for every 1 < i < d\ 
{ dxi J 



and for u,v G J 7 , 



A(d,a)a a , , 
dxdy. 



£(u,v)=[ Vu(x) ■ Vv(x) dx + ^ / (u(x) - u(y))(v(x) - v(y)) 

In probability theory, the quadratic form {£, W l,2 (R d )) is called the Dirichlet form of X a . A 
statement is said to hold quasi-everywhere (q.e. in abbreviation) if there is a set N having zero 
capacity with respect to (£%, W 1,2 (M. d )) such that the statement holds everywhere outside N. Here 
£i(u, u) := £{u,u) + j Rd u(x) 2 dx. The function J a (x,y) := a a A(d, a)\x — ?/| _ ( d + Q ) is the Levy 
intensity of X a . It determines the Levy system for X a , which describes the jumps of the process 
X a : for any non- negative measurable function / on M + x M. d x R d , x G M. d and stopping time T 
(with respect to the filtration of X a ), 



f(s, Xf_,X%) 

s<T 



f(s,Xf,y)J a (Xf,y)dy)ds 



(1.4) 



(see, for example, [181 Proof of Lemma 4.7] and [HI Appendix A].) 

The purpose of this paper is to establish the scale invariant version of the boundary Harnack 
principle in Theorem 11.41 To state this theorem, we first recall that an open set D in M. d (when 
d > 2) is said to be C 1 ' 1 if there exist a localization radius R > and a constant A > such that 
for every Q G dD, there exist a C^-function <fi = 4>q '■ M^" 1 — > R satisfying 0(0) = V0(O) = 0, 
||V0||oo < A, |V(/>(x) — V0(y)| < A|a; — y\, and an orthonormal coordinate system CSq: y = 
(yi, ■ ■ ■ i Vd-uVd) ='■ (y, yd) with its origin at Q such that 

B(Q, R) n D = {y = (y, y d ) G B(0, R) in CS Q : y d > cj>{y)}. 

The pair (R, A) is called the characteristics of the C > open set D. By a C ' open set in R we 
mean an open set which can be written as the union of disjoint intervals so that the minimum of the 
lengths of all these intervals is positive and the minimum of the distances between these intervals 
is positive. Note that a C ' open set can be unbounded and disconnected. 
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For any x E D, let 5d{x) denote the distance between x and dD. It is well known that any 
C > open set D satisfies the uniform interior ball condition: there exists R < R such that for every 
x E D with 5d(x) < R, there is Q x E dD so that \x — Q x \ = 5d(x) and that B(x,R) C D for 
x = Qx + R(x — Qx)/\ x — Qx\- Without loss of generality, throughout this paper, we assume that 
the characteristics (R, A) of a C ' open set satisfies R = R < 1 and A > 1. 

For any open set D C M d , := inf{t > : X® ^ D} denotes the first exit time from D by X a . 

Definition 1.1 A real- valued function u defined on R d is said to be harmonic in D C R d with 
respect to X a if for every open set B whose closure is a compact subset of D, 



E, 



u(X? a )\ < oo and u(x) = E x u(X° 



for q.e. x E B. 



(1.5) 



Note that by using the Levy system of X a , we have 



E, 



> E 



E 




u(y) I , ^ — rn — ] (is 

B Vy;i \X%-y\ d +« J 



Hence if n is a harmonic function in D with respect to A a , then A integrable 

on 5 C for any relatively compact open subset B with B C D. It follows from Theorems 1.2 
and 1.3 of [20] that all harmonic functions in D with respect to A a are continuous on D, since 
every harmonic function in D with respect to X a can be approximated locally uniformly in D by 
functions that are bounded on M d and harmonic with respect to X a in relatively compact open 
subsets of D. Therefore, for any harmonic function u in D, (| 1 . 5 j) holds for every point x E D. The 
above also implies that any harmonic function u in D with respect to X a is locally bounded in D 
with J Rd \u(y)\(l A \y\~^ d+a ^)dy < oo. A function u is said to be in W^(D) if for every relatively 
compact subset B with B C D, there is a function / E VF 1,2 (M d ) such that u = f a.e. on £?. The 
following analytic characterization of a function u being harmonic in D with respect to X a follows 
immediately from Example 2.14 in |16j . 



Proposition 1.2 Let D be an open subset ofM. d . Then the following are equivalent. 

(i) u is harmonic in D with respect to X a ; 

(ii) u is locally bounded in D, f Rd \u(y)\(l/\\y\~( d+a ^)dy < oo, u E W^(D) and (A+a a A a / 2 )u = 

in D in the distributional sense: for every <f> E C£°(D) 



V«(x) • Vcp(x) dx + 



(u(x) - u{y))(4>{x) - 4>(.y)) ^ d,a } d + a dxd v = °- 



f - y\ 



The following uniform Harnack principle will be used to prove the main result of this paper. 

Proposition 1.3 (Harnack principle) Suppose that M > 0. There exists a constant Cq = 
Co(a,M) > such that for any r E (0,1], a E [0, M], xq E M. d and any function u which is 
nonnegative in M. d and harmonic in B(xo,r) with respect to X a we have 

u(x) < Cou(y) for all x,y E B(xo,r/2). 
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Let Q G dD. We will say that a function u : M. d — > R vanishes continuously on D c n -B(Q, r) if 
n = on D c n B(Q,r) and n is continuous at every point of dD n B(Q,r). The following is the 
main result of this paper. 

Theorem 1.4 Suppose that M > 0. For am/ C 1 ' 1 open set -D inR rf uni/i the characteristics (R,A), 
there exists a positive constant C = C (a, d, A, R, M) such that for a G [0, M], r G (0, R], Q G <9.D 
and any nonnegative function u in R d t/iaf is harmonic in D n B{Q,r) with respect to X a and 
vanishes continuously on D c n -B((J, r), we /iawe 

^4<C^44 for every x,yeDr\B{Q,r/2). (1.6) 

When a changes from to A/, A + a°A a / 2 changes continuously from A to A + M a A Q / 2 . So 
Theorem 11.41 savs that the BHP holds uniformly for the family {A + a Q A Q//2 , a G [0, M]} of pseudo 
differential operators in the sense that the constant C in (jl.6p can be chosen to be independent of 
a G [0, M]. Note that a = corresponds to the classical case of the boundary Harnack principle 
for the Laplacian. We will therefore in the rest of the paper assume that a G (0, M]. 

As far as we know, this is the first time that a BHP has been established for non-local integro- 
differential operators that have second order differential operator components in dimension two and 
higher. Unlike (11. ip and the paragraph following it, in this paper we are concerned with the above 
BHP for C > open sets only. The main focus and goal of this paper is to get the explicit decay rate 
of harmonic functions near the boundary of D as in (|1.6p and to show that the BHP is uniform in 
a G [0, M\. We emphasize that (jl.6p is not true in Lipschitz domains even in the classical case of 
BHP for the Laplacian. However, a uniform Carleson type estimate is shown to hold for Lipschitz 
open sets in Theorem 14.31 The BHP of above type is very useful in studying other fine properties 
of the process. For example, we will use it to derive sharp Green function estimates of X a in C 1 ' 
open sets in a forthcoming paper [17] . 

For a > 0, X a and X := X 1 are in fact related by a scaling. More precisely, for a G (0, M], 
X a has the same distribution as AX^-2 t , where A = a a ^ a ~ 2 ^ > M a ^ a ~ 2 \ Consequently, if u is 
harmonic in an open set U with respect to X a , then v{x) := u(Xx) is harmonic in \~ l U with 
respect to X. Hence the uniform Harnack inequality of Proposition 11.31 follows from the Harnack 
inequality for X. The latter is known, see Theorem 6.7 of [20] or Theorem 4.5 of [38] . However 
the uniform BHP of Theorem 11.41 can not be obtained by such a scaling argument from the BHP 
of X. This is because for a C ' open set D with the characteristics (R, A), X~ 1 D is, in general, a 
C ' open set with C ' characteristics (R/X, AA), which tends to (0, oo) as A — ► oo. 

For each fixed «o G (0,2), when a changes from ao to 2, the operator A + a°A a / 2 evolves 
continuously from A + a a °A a °/ 2 to (l + a 2 )A. So in view of Theorem ll.4] it is reasonable to expect 
that one can get the BHP for A + a a A a ^ 2 uniformly both in a G (0,M] and in a G [qo,2). We 
believe this is the case and that it can be achieved by carefully keeping track of all the comparison 
constants in the arguments of this paper. However in order to keep our exposition as transparent 
as possible, we are content with establishing the results stated in Theorem 11.41 and leave the details 
of the proof for the last claim to interested readers. 

Our method of establishing the above BHP is different from those in UD] for symmetric 
stable processes and in [31] for more general subordinate Brownian motions. The reason that the 
approaches in |4"0~1 [3~T] do not work well in our setting lies exactly with the fact that X a leaves 
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open set D by jumping out across the boundary dD as well as by continuously exiting D through 
the boundary of D. To circumvent this difficulty, in this paper we adopt the ideas from [8] for the 
BHP of censored stable processes, which are further refined in [26] . That is, we use suitably chosen 
subharmonic and superharmonic functions of the process X a (or equivalently, of A + a Q A a / 2 ) to 
derive some exit distribution estimates that are needed to establish the BHP. However, had we 
done it in this way directly, we would only get the BHP for A + o Q A Q//2 with a G (1,2). The 
reason is that, when D = W\_ := {x = (an,-- - ,x d ) G R d : x 1 > 0}, we need to consider testing 
functions w p (x) = {x\ V 0) p for p > 1. But for w p to be A a//2 -differentiable in H+, see ljl.3|) . one 
requires p < a, which would be impossible when a G (0, 1]. To overcome this difficulty, for each 
A > 0, we consider the finite range (or truncated) symmetric a-stable process Y A obtained from Y 
by suppressing all its jumps of size larger than A. The infinitesimal generator of Y x is 

■= lim / («(y) - <*))t^L d y- ( L7 ) 

e l° J{y<=R d :e<\y-x\<\} \ x ~ ^ 



When A = 1, we will simply denote A^ 2 by A^ 2 . Then w p is A^ 2 -differentiable in for 
every p > 0. Observe that X a := X° + aY l l a is a Levy process obtained from X a = X° + aY 
by suppressing all its jumps of size larger than 1 and that the infinitesimal generator of X a is 
A + a Q A^ //2 . From this, we can obtain suitable exit distribution estimates for the Levy process X a . 
The desired estimates for X a can then be obtained from that for X a by adding back those jumps 
of X a of size larger than 1. Such an idea has already been used in [21] to study Schramm-Lowner 
evolutions driven by one-dimensional symmetric stable processes. We remark that the BHP in 
Theorem 11.41 for the case of a = 1 has also been mentioned in Remark 5.2 of Guan [26]. However, 
no precise statement (such as the range of a) nor a proof is given in that paper. 

The rest of the paper is organized as follows. In Section 2, we derive estimates on A^ 2 w p . These 
estimates are then used in Section 3 to obtain exit distribution (or harmonic measure) estimates for 
the finite range process X a and then for the desired process X a . In Section 4, we first give the proof 
of Proposition 11.31 and then establish a Carleson estimate for non-negative harmonic functions of 
A + a a A a / 2 in Lipschitz open sets. Then using these results, the proof of Theorem 1 1.41 is presented. 

Throughout this paper, we use the capital letters C\ , C2, • • • to denote constants in the statement 
of the results, and their labeling will be fixed. The lowercase constants ci,P2, •■■ will denote 
generic constants used in the proofs, whose exact values are not important and can change from 
one appearance to another. The labeling of the constants c±,C2,--- starts anew in every proof. 
The dependence of the constant c on the dimension d > 1 and a G (0, 2) may not be mentioned 
explicitly. The constant M > will be fixed throughout this paper. We will use ":=" to denote a 
definition, which is read as "is defined to be". For a, b G R, aAb := min{a, b} and a V6 := max{a, b}. 
For every function /, let / + := / V 0. We will use d to denote a cemetery point and for every 
function /, we extend its definition to d by setting f(d) = 0. We will use dx or m^dx) to denote 
the Lebesgue measure in R rf . For a Borel set A C R d , we also use \A\ to denote its Lebesgue 
measure and diam(^4) to denote the diameter of the set A. 
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2 Truncated fractional Laplacian estimates for power functions 

In this section, we give some estimates which will be used later. Recall that the fractional Lapla- 
cian A°/ 2 and the truncated fractional Laplacian A^ 2 := A^ 2 are defined in (|1.3p and (jl.7p . 
respectively. 

Lemma 2.1 For x 6 M. d and p > 0, set w p (x) := (x^) p . Then there are constants R* £ (0, 1), 
C\ > C2 > depending only on p, d and a such that for every x£K' i with x\ £ (0, i?*] 

|A^ /2 Wp (x)| < d /orp>a, (2.1) 

\L a J 2 w v {x)\ <Ci|logxi| /orp = a, (2.2) 

C 2 x\- a < A a J 2 w p {x) < C x x\~ a for a/2<p< a, (2.3) 

- Ci < A^ /2 u; p (x) < -C 2 /or p = a/2, (2.4) 

and 

- Cix?" < Ad /2 ^ P (^) < -C 2 x{~ a for0<p< a/2. (2.5) 

Proof. First note that using integration by parts and a change of variable, we get that for p, x > 
and g G (0, l/(x + 1)), 

l-e p _ 1 1 /-1-e 

dz = - I (z p - l)d(l - z)" a (2.6) 



(1 - z)^ 1 a J 



1 



a 



Q: 



(1 - e) p - 1 1 p Z" 1 ^ zP" 1 ^ 



and 



ae" a a J (1 — z) 



dz = -- / (z p - l)d(z - l)~ a (2.7) 



1+£ (z-l)^ 1 ay 1+ . 



1 



a 



31+1 



1+e a Ji + 



x+1 



(z p - l)(z - l)~ a x +- I z p -\z-l)- a dz 



-e 



= ii±£l! — - + _ -( x + + £ ^-p-^i - z)- a dz. 

ae a a a a J * 

x + 1 

For p > and x G (0, 1), by a change of variable 

a«/2 / \ An \ v f w p(v) ~ w p( x ) A 
A/ w p (x)=A{l,a) hm J^—^—^—i {e<y _ x ^ 1} dy 

rx+l y P_ x P r0 dy 



rx+l yP _ X P r 

A(l, a) hm J \ x _ y \i+a Hv-xfr*} d y - A (h a)x p J 



x - y\ 1+a 



x-1 





x + 1 

A(l,a)x p - a \imJ^ Z _ ~l +a l { | z _i| >£M dz - A(l, a)x p J {x-yY^dy 



i_ £ „ . .«±i 



^(l,«)^lim| / * d* + [ X , Z \* dz)- A ^(x p - a -x p 



eio U J 1+e (z / a 
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So we have by (|2l)]) - (|277|) that for p > and x G (0, 1), 



(2.8) 



4(_^) lim U + (i- e y + (i + e y-2\ _ 41^ {xP _ a + 4(1^ 



.4(1, a) 



a 



(x + 1)* 



4(1, a)p 



-x 



p-a 



O 



1 z a- p -l _ z p-l f-^j z p-l 

-dz 



4(1, a) 



a 



2x p - (x+l) p +px p " 



X + l 



1 ^Q— p— 1 _ z p— 1 

(1-^) Q 

.1' + 1 



o (l-*) a 



dy 



dz 



x+l z- 



p-i 



o (l-^) a 



dz 



Note that for p > a, 



sup x 

xe(o,i 



p-a 



1 ^.a— p— 1 _ ^.p— 1 



(1-zy 



dz 



x+l z- 



p-i 



o (1 " z) 1 



-dz 



< oo. 



So when p > a, 



sup I A^u^x)! < oo. 

x6(0,l) 



When p = a, there exists an r, > such that for < x < r* 



(2.9) 



1 z" 1 - z a ~ x 



dz 



(1 " A a Jo (1 - *) a 



x+l Z 



a-1 



-dz 



< 



1 z- 1 



(1 - z)' 



-dz< (l + r*) a log((l + r*)/x). (2.10) 



It is easy to see that 



sup 

x6[r«,l] 



1 -1 _ a-1 



(1 " z) a Jo (1 " z) 



x + l Z 



a-1 



-dz 



< OO. 



(2.11) 



On the other hand, when p G (0, a), 



sup 

xe(o,i] 



1 z a— p— 1 1 



(1 " *)< 



x+l z 



p-i 



o (1 " z) a 



-dz 



< oo. 



(2.12) 



As 



lim 

x^0+ 



1 y Ot— P— 1 _ 1 



X+l 



(l-z)< 



while for p = a/2, 



1 ^.a— p— 1 _ Z P—1 



dz 



*+i z 



p-i 



(l-z)<* 



dz 



> if p G (a/2, a) 

<0 ifpG(0,a/2) 



x+l z 



p-i 



-dz 



o (1-^) Q Jo (l-«) 



-dz, 



we conclude from (|2.8j) - (12.12|) that there are constants r\ G (0, 1) and C\ > C2 > depending on 
p and a so that when p = a, 



| A^u^x) < Ci|logx| for x G (0, ri] and sup | A° /Z w p (x)| < 00, (2.13) 

xG(ri,l) 



-a/2 
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when p = a/2, 



d < A" 7 w p {x) < -C 2 for x G (0,n] and sup |A° 7 w p (x)\ < oo, (2.14) 

a;6(ri,l) 



when p G (a/2, a), 

C 2 x p ~ a < A" /2 w p (x) < Cix p ~ a for x G (0,n] and sup |A" /2 w p (x)| < oo, (2.15) 

x€(rt,l) 

and for p G (0, a/2), 



Cix p ~ a < A" /2 w p (x) < -C 2 x p ~ a for x G (0, n] and sup |A" /2 u; p (x)| < oo. (2.16) 

xe(ri,l) 

On the other hand, for x > 1, 



Ai /2 w p (x) 



/•— yP - 1 
Vl(l,a)x p "Ijmy^ | y _!|i+a 1 {|y-i|> e }^ 



1 ~ e yP-1 , yP-1 



^(l,.)^f ^ 1 + ^^-^- 2 ^. 
Note the above integrand 

(l + n)P + (l-n)P-2 
u 1+a 

is of the order u l ~ a near zero. So for p > and a G (0, 2), there is a constant Co = co(p, a) > so 
that 

|A" /2 w p (2;)| < c x p " 2 for x > 1. (2.17) 

With ri G (0, 1) as in (|2.13p - (|2.16p . the above inequality in fact holds for x > r\. 

The estimates (|2.9p - (|2.16p prove the Lemma in dimension d = 1. Now we consider the case 
d > 2. For each fixed x G we use the spherical coordinates 

(j/i, ... , yd) := x + (r cos #1, r sin6>i cos 6 2 , . . . ,r sin#i . . . cos 0<j_i, r sin#i . . . sin#<2_i) 

where r > 0, < 0\, . . . , 9d- 2 < vr and < 0d-i < 2tt. Let 

<P(9) : = 0(0 l5 . . . , 6 d _ 2 ) := sin^ 2 9 1 sm d ~ 3 9 2 ... sin d _ 2 . 

Then for x G M. d with xi > we have 

lim/ ((2/i + ) P -x?), 



|x — y\ d + c 



2?r 



tt/2 riT rn r-2n ^ rl (( r _|_ 



aim j dfl! y ^ 2 --- y o de d . 2 ( f>(e)do d -i(c OS e 1 f j - cosdl rl ; a cosdl ' dr 

fir fir f-n r-2n ft (( r x l ) + )P ( x l )P 

+ iim dej de 2 -- de d _ 2 <j ) {e)de d _ 1 {- C ose l y ^ cos(>1 \L cosdl dr. 

40 7^/2 JO Jo Jo Je r + 
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By the change of variable r = t — x\/ cos9\ for 9 G [0,7r/2) and r = —t — x\/ cos9\ = — t + 
xi/ cos(7r — #1) for 9 G (7r/2,7r], we get 

lim / ((y+f - x\) 



£ 



TO J{ym d :l>\y-x\>e} yyal ' U \x ~ y\ d+a 

rn/2 i"K fir fits rl-\ — (t + ) p — ( Xl ) v 

hm / d9 1 d9 2 -- d9 d _ 2 / ^^(cosfli)* / C ° S 1 , ; . C °l e \' dt 

^ Jo JO Jo Jo Je+-^r- * ~ Xl ' COS 01 1 

cos 

2 - - '-«+=l2*7 (i + ) p -(^Ty) p 



+ hm / d9 1 d9 2 --- d6 d - 2 / mddd-i COS 7T-0! f / , / ll+a ^ 

Jo Jo Jo J-i+ , X1 „ , I* — a?i/cos(7r — 6' 1 )| 1 +« 



COs("7T — ) 



hm / d0! d9 2 -- d9 d „ 2 / ^de dr . 1 (cose 1 y / cos 1 =f . „lL a dt 

dOJo Jo Jo Jo Je+^s- \t-Xl COs9 1 \ L+a 

' COS fl^ 

+ hm / d9 1 d9 2 ..J d9 d „ 2 / 0(0)^ (cob / ^ 1 . V . ™ 6l ' dt 
do Jo Jo Jo Jo J-i+^V I* — ^l/cos 6> 1 | 1 +« 

COS f]^ 

hm / (t+)P - ( - 

e|0 7 {teR:1> | t __£L_| >£} |t- Xxl COS I 



/•TT/2 /-7T f-TT [2-K I r- (t + ) P — ( Xl ) P \ 

Jo Jo Jo Jo V ei0y {t6R:1> u _fi| >e> |t- X 1 /cOS^i| 1 + a / 



Therefore we have 
A^ /2 w p (a;) 



a) r de 2 • • • r d9 d ^ 2 r ^d^a* e^A^wp ( xi 



-4(1, a) Jo Jo Jo Jo \cos9x 

Aid a\ / , arccos(ii/ri) /-tt /-tt /-27T _ / \ 

/ <#i / d6 2 --- / dfl d _ 2 / WW^coBOipA^Wp 

"4(1, a) Jo Jo Jo Jo V cos6l i/ 

+4tT 1 4 ^ /"W" f d0 d _ 2 /^^^(cos^A^f^-) 

-4(1, a) Jarccos^i/n) Jo Jo Jo \cos9xJ 



The conclusion (|2.ip - ()2.5p now follow immediately from the above equality and the estimates (|2.9p - 
(|2.17p . where we use (|2.17p to bound the second integral above by c\x\/r\ for some positive constant 
ex. □ 

Remark 2.2 A careful evaluation of (|2.8p in fact shows that lim x ^o+ A^ 2 w p (x) = c / when 
p > a. At the hrst glance, this may look surprising, as in the Brownian motion case (which 
corresponds to a = 2), Ax^ = p(p — l)xj'~ 2 . The bound in (|2.ip is due to the non-local nature of 
the operator A^ 2 for a G (0, 2). However a more careful analysis of (|2.8p reveals that for p > 0, 

A° /2 w p (x) x (2- a) ((p-a)' 1 - l) + p(p - l)x p ~ a for x G (0,n) 

as a | 2. It is not difficult to see that as a f 2, A" //2 w p (x) converges to Aw p (x). □ 

Recall that for A > 0, the operator A^ 2 is defined by (jl.7p . Note that 

A^(x) = A- a (A^(A-))(A" 1 x). (2.18) 
Thus, from Lemma |2. II and (|2.18p . we get the following corollary. 
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Corollary 2.3 For x G M. d and p > 0, set w p (x) := (xf) p . Then there are constants -R* G (0, 1/2), 
Ci > C-2 > depending only onp, d and a such that for every A > and x G K d with x\ G (0, \R*), 

\K a £w p {x)\ < dX^ for p>a, (2.19) 

\A a d ^w p (x)\ < dllog^x/A)!, /or p = a, (2.20) 

C 2 x?~ a < A^w p (x) < dx?-* for a/2<p< a, (2.21) 

- CiA- Q / 2 < A^w p (x) < -C 2 \- a ' 2 for p = a/2, (2.22) 

and 

- C lX p ~ a < A^w p (x) < -C 2 x p ~ a for0<p< a/2. (2.23) 

3 Estimates on harmonic measures 

Recall that for any open set P C I'', = inf{i > : Xf G' U} is the first exit time from U by 
X a . 

Lemma 3.1 For every b G (0,oo), there exist C3 = C^{M,b) > and C4 = C±(M, b) > such 
that for every x G R d , a G (0, M] and r G (0, 6], 

C 3 r 2 < E X0 [r^ xo>r) ] < C 4 r 2 . (3.1) 

Proof. See Lemmas 2.2 and 2.3 in |37j or Lemmas 2.3 and 2.4 in [20J for a proof. □ 

In the remainder of this section, we assume D is a C ' open set with characteristics (R, A). 
Recall that we are always assuming that R < 1 and A > 1. For notational convenience, throughout 
the rest of this section, we put 

R 



4Vl + A 2 ' 
Define 

p Q (x) := x d - 4>q(x), 

where (x, xj) is the coordinates of x in CSq. Note that for every Q G <9-D and x G -B(Q, R)n D we 
have 

(l + A 2 )- 1 /^^) < J D ( X ) < PQ ( X ). (3.2) 
Recall that R* is the constant in Lemma l2.1i 

Lemma 3.2 Fix Q G dD and the coordinate system CSq so that 

B(Q, R)nD = {y= (y, y d ) G B(0, R) m CSq : y d > <j>{y)} . 

For p > a/2, let 

h p(y) ■= (pQ(y)) p iDnB(QAr )(y)- 

Then there exist Ci = Ci(a,p, A, R) > Q, i = 5, 6, 7, independent of the choice of the point Q G dD 
such that 
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(i) in the case ^ < p < a, for all x £ D such that pq(x) < ro A R* and \x\ < tq, we have 

C 6 { PQ (x)f- a < A« /2 h p (x) < C 5 ( PQ (x)f- a ; (3.3) 

(ii) in the case p > a, for all x 6 D such that pq(x) < ro A R* and \x\ < ro, we have 

\A a d /2 h p (x)\<C 7 ; (3.4) 

(iii) in the case p = a, for all x £ D such that pq(x) < ro A R* and \x\ < ro, we have 

\A a d /2 h p (x)\<C 7 \log(p Q (x))\. (3.5) 

Proof. In this proof our coordinate system is always CSq. Fix x = (x,x d ) £ D such that 
Pq{x) < ro A i?* and \x\ < ro, and choose a point xq £ dD satisfying x = Xq. Denote by 
li(xo) the inward unit normal vector at xq for dD and set Q(y) = (y — xq,Tl(xq)} for y £ R d . 
Then IL := {y : &(y) = 0} is the hyperplane tangent to dD at the point xq. The function 
r* : R rf_1 — > R describing the plane II is given by T*(y) = c/)q(xo) + V^q(xo)(2/ — ^o)> and it holds 
that ( (y,T*(y)) - x , "ra (z )) = 0. We also let 

A := {y : T*(y) < y d < 4> Q (y) and \y - x\ < r } |J {y : T*(y) > y d > </> Q (y) and \y - x\ < r }, 
E := {y £ D\ A : \y-x\ < r and p Q (y) < r (2 + A)}. 

Note that, if \x — y\ < ro and y £ D, 

pq{y) < \yd - xd\ + \xd - 4>q(x)\ + \<t>Q{y) - 4>q(x)\ < r {2 + A). 
On the other hand if \y — x\ < ro and PQ(y) < ro(2 + A), then 

\y\ 2 = \V\ 2 + \y d \ 2 < (2r ) 2 + (r (2 + A) + \4> Q (y)\) 2 < ^ttf) * < ^ 

Consequently, we have 

Dr\B(x,r ) cDn{y: \y-x\<r and p Q (y) < r (2 + A)} C DC\B(0,R). (3.6) 

Let h(y) := h x (y) := (y d — T*(y)) + for y £ M. d . Since V(/)q(x) = VP*(af), by the mean value 
theorem and the C ' condition on c/)q, 

\h(y)-P Q (y)\ = \<j> Q (y)-r*(y)\ (3.7) 
= \<I>q(v) -4>q(x) - V0 Q (i) • (y-x)\ <A\y-x\ 2 , y £ E. 

For y £ R d , define 5 n (y) := dist(y, n) and D r * = {y £ R d : y d > T*(y)}. Let 

b x :={l + \V^ Q (x)\ 2 ) 1/2 and h x , p (y) := (h{y)f for p > a/2. 



Note that 1 < b x < \/l + A 2 and h XtP (x) = h p (x). 
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Recall that R* and C\ > C2 > are the constants in Lemma [2. li Since h(y) = b x 5 n (y) on D^*, 
by Lemma \2.1\ it holds that for y £ L>r* and 5 n (y) < i?*, 

C2 KK(y)) P ~ a < K' 2 h x ,viy) = KK'\^{y)Y < d Wn(l/)) p - a when a/2 < p < a, (3.8) 



K /2 K P (y)\ = K\K'\^{y)) P \ < Ci < Ci(l + A 2 f/ 2 when p > a, (3.9) 



-a/2 



|Af^, P (y)| = ^.|A^ 2 (5 n (y)n < Ci^llog^d/))! 

< Ci(l + A 2 ) p / 2 |log(5 n (y))| whenp = a. 



(3.10) 



Note that &a;<5 n (x) = Pq(x). Applying (|3.2p and (|3.8p to the point x gives that, for a/2 < p < a 



C2 PQ {xy- a < c 2 bz PQ { x y- a < A a d /2 h x , p (x) < c^pQixy-* < d(i + at /2 pq( 



Note that by (1331) . 



(3.11) 



(3.12) 



lim 

E. 



= «4(cZ, a 
< .A(d,a 

+^4(d, a) lim 



(/i p (y) - h p>x (yj) 



^ J{l>\y-x\>e} \x-y\ d+a 
(h p {y) - h PtX {y)) 



dy 



dy 



{i>\y-x\>r } \x-y\ d+a 

\h P (y) - h P:X (y)\ 



dy 



£° J{r >\y-x\>e} \x - y\ d+c 

(h p {y) - h p , x (y)) 



dy 



{i>\y-x\>r } \ x ~y\ 



d+a 



\h P (y) - h PtX (y)\ 



= : h+h + h. 
We claim that, if p > a/2, then 



h + h + h< c 



(3.13) 



for some constant cq = co(a,p, A, R). Together with (|3.9p - (|3.12p this will establish the desired 
estimates (|3,3p - (|3,5p with constants depending on a, p, A and R. 

Clearly I± is bounded by some positive constant. 

For y G A, we have 



\h x ,p(y)\ + \h P (y)\ < \vd-r*(y)F + \vd-<f>Qm p < 2|0o(y)-r 

< 2|0o(5)-0o(x)-V0o(x)-(i/-5)|f < 2AP|y-5?| 2 P. 
Furthermore, since, on {|y — x\ = r < ro}, \<Pq(jj) —T*(y)\ < A\y — x\ 2 = Ar 2 , 

rnd-i ({y ■ \y-x\ = r,T*(y) < y d < <j) Q (y) or T*(y) > y d > <t>Q(y)}) < c\r a 



(3.14) 
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for some constant c\ > if r < tq. This together with (|3.14p yields that 

t <r ma \ f° f i ( >*, P (s/)l + IMs/)l (j w 

h < A{d,a) l A (y) — \~ _~ H+a m d ^ 1 (dy)dr 

JO J\y-x\=r W x \ 

rro 

< c 2 r~ d+2p ~ a m d _i({y £ A: \y — x\ = r})dr 

Jo 

rro 

< cic 2 / r 2p ~ a dr < c 3 . 

Jo 

Note that for y £ E 

My) - h XjP (y)\ =\(h(y)T - ( PQ (y)) p \ < c 4 (h(y))^-\h(y) - p Q (y)\, (3.15) 
where (p — 1)_ := (p — 1) A 0. In the last inequality above, we have used the inequalities 

\y°-a p \ < &P -1 |&-a| for a,b > 0, <p < 1 

and 

\V ~ a p \ <(p+l)\b- a\ for a, b 6 (0, 1), p > 1. 

For y = (y, y d ) £ K d , we use an affine coordinate system z = (z, z d ) to represent it so that 
z d = Ud~ r*(y) and z are the coordinates in an orthogonal coordinate system centered at xq for 
the (d — l)-dimensional hyperplane II for the point (y, T*(y)). Denote such an affine transformation 
y \—* z by z = 9(y)- It is clear that there is a constant C5 = cs(A, R) > 1 so that for every y £ M. d , 

<^ l \y ~x\ < \z\ < c 5 \y-x\, c^\y - x\ < - V(x)\ < c 5 \y - x\ 

and that 

9(E) C{z= (z, z d ) £ R d : \z\ < c 5 r and < z d < c 5 r }. 

Denote x d — F*(x) by w; that is, 9(x) = (0,w). Hence by (|3.7p and (|3.15p and applying the 
transform 9, we have by using polar coordinates for z on the hyperplane II, 

%) r"lr^w ,r -'" |Sf * 

E \y~x\ d + a U ~ J nE) \ z - (p, W )\d+a 



csro , / pc 5 ro r d-2 

(P-I)- 



- C8 (, Zd " [Jo (r + \z d - W \)^ dr I dZd 



csro <p-l)_ / 1 1 



~ ° 9 Jo Zd Vl^-^l"" 1 (c 5 r + \z d -w\)°<- 1 ) dzd 

< cio / n -vt < en < 00, 

■A) 4 u;|a-l 

where all constants depend on a, p, A and i?. The last inequality is due to the fact that since p > 0, 
< a < 2 and (1 — p) + + a — 1 = max{a — p, a — 1} < 1, by the dominated convergence theorem, 

rCBro I „ , „„|l-a. 



(w 



) := Jq 5T ° z d \z d — w\ l a dz d is a strictly positive continuous function in x d £ [0, c^tq] and 
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hence is bounded. Thus we have proved the claim (|3.13p . hence completing the proof of the lemma. 
□ 

Since D is a C > open set with characteristics (R,A), for every A > 1, XD is a C ' open set 
with uniform characteristics (R, A). Thus, by the previous lemma and (|2.18p . we get the following 
as a corollary. 

Corollary 3.3 Fix Q G dD and the coordinate system CSq so that 

B(Q, R)C\D={y= (y, y d ) G B(0, R) in CSq : y d > <f>(y)}. 

Let 

Mv) : = (pQ{y)) Pl DnB(QAr )(y)- 

Then there exist Ci = Ci(a,p, A, R) > 0, i = 5, 6, 7, independent of the choice of the point Q G dD 
and A > 1 such that 

(i) in the case ^ < p < a, for all x G D such that pq(x) < ro A i?* and \x\ < ro, we have 

C 6 { PQ (x)f- a < A^hp(x) < C 5 ( PQ (x)f- a ; (3.16) 

(ii) in the case p > a, for all x G D such that Pq(x) < ro A i?* and \x\ < ro, we have 

-a/2, , 

(iii) in the case p = a, for all x G D such that pq(x) < ro A R* and \x\ < ro, we have 

\A a d ^h p (x)\<C 7 \log(p Q (x)/X)\ . (3.18) 

The following scaling property of X a will be used below: If {X^' D , t > 0) is the subprocess in D 
of the independent sum of a Brownian motion and a symmetric a-stable process on M rf with weight 
a, then (AX?!^.,t > 0) is the subprocess in XD of the independent sum of a Brownian motion and 
a symmetric a-stable process on M d with weight aX^ a ~ 2 ^ a . So for any A > 0, we have 

V\d ~ 2)/a (*, x, y) = X~ d p a D (X- 2 t, X^x, X^y) for t > and x,y £ XD. (3.19) 

By integrating the above equation with respect to t, we get 

G a x x ^ 2)/a (x,y) = X 2 ~ d G a D (X~ l x, X~ x y) for x,y e XD (3.20) 

where 

/•oo 

G a D (x,y) := / p a D (t,x,y)dt 



\A a d 'x%(*)\ < C 7 X p ' a ; (3.17) 



is the Green function of X a in D. It is well known that the Levy measure of X 1 has the intensity 

J\x,y) = j\\x - y\) = A(d, a)\x - y\^ d+a \ 
Thus by a scaling argument, we get that the Levy intensity of X a is 

J a (x,y) = f(\x - y\) = a a A(d, a)\x - y\^ d+a \ 
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which gives the Levy system (jl.4p of X a . 

By a A-truncated symmetric a-stable process in M. d we mean a pure jump symmetric Levy 
process Y x = (Y t x ,t > 0,P X ) in R d with Levy density a)\x\- d - a 1{| X |<A}- Note that the L «vy 
exponent ip x of y A defined by 



is given by 



^\i)=A{d,a) 



for every x G M d and £ G M d , 
1 - cos(£ • y) 



\d+a 



-dy. 



(3.21) 



>{\ y \<\} \y\ 

Suppose that yV a is a (A/a)-truncated symmetric a-stable process in M. d which is independent 
of the Brownian motion X°. For any a > 0, we define 



X 



a, A 



X? + a ? t X/a , t>0. 



Note that from (|3.2ip we can easily check that for any b > 0, 

i) x (b£) = b a ^ Xh (^) for every £ G R d . 

Thus for any a > and 



(3.22) 



E, 



-*l€l a E, 



-t(|£|W»(«0) 



-t(|£| 2 +a»V A (0) 



Therefore JT a ' A has the same distribution as the Levy process obtained from X a by removing jumps 
of size larger than A. The above observation also gives us that the infinitesimal generator of X a,x 
is A + a a A^'?, and the Levy intensity for X a ' X is 

J a > X (x,y) := a a A(d,a)\x - y\~^ l {lx - yl<x} ■ 



The Levy intensity describes the jumps of the process X a,x through the Levy system: for any 
:abl< 

to the filtration of X a > x ) 



non-negative measurable function / on M + x M. d x R rf , x G M. d and stopping time T (with respect 



E, 



j2f(s,x:?,x^ 



s<T 



E. r 



f(s,X«> x ,y)J a ' X (X^ x ,y)dy)ds 



(3.23) 



For any open set U C M. d , let t^' = inf{t > : X^' ^ U} be the first exit time from U by X a ' X , 
and denote by X a ' X,u the subprocess of X a,x killed upon leaving U. When A = 1, we simply write 
X a for J"' 1 and for t^ 1 . The following scaling property will be used in the next lemma: by 
(|3.22p . we see that for every A, a, b > and £, x G M. d , 



E, 



D iHKx a b L\ t -x a ' x )) 



A/a 



b-2« r > 



Thus, if {X^ ,X,D ,t > 0} is the subprocess of {X^ ,A ,t > 0} in D, then {bX^^ ,t > 0} is the 
subprocess of {X" 6 ' > 0} in bD. In particular, if > 0} is the subprocess of 



e -m 2 e -b- 2 t^/-(abO = e -f(| ? |2 +a -fe«-2^ A (0). 



>a,X,D 
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{X?,t > 0} in D, then {XX^ 2f ,t > 0} is the subprocess of {X^ 2)/a ' X ,t > 0} in XD. So for any 
A > 0, we have 

p^ (a ' 2)/a ' X {t,x,y) = A^'V^t.A'^.A- 1 !/) for t > and x,y G AD (3.24) 

where p^ ,X (t,x,y) is the transition density of X a ' A,D . By integrating the above equation with 
respect to t, we get 

G^ a ~ 2)/a ' X (x,y) = \ 2 - d G a /(\- l x,\- l y) for x,y G XD (3.25) 

where 

f oo 



G~D X (x,y) := / p£' X (t,x,y)dt 
Jo 



is the Green function of AT a ' A in D. 

For our reader's convenience, we summarize some notations below. 



Process 


Generator 


Levy (jumping) kernel 


X° 


A 





Y 


A a/2 


A{d,a)\z\- d - a 


aY 


a a A a/2 


a a A(d,a) \\z\~ d - a 


Y x 




A(d, a )\zr d ~ a l {lzl<x} 


X a := X° + aY 


A + o a A"/ 2 


a a A{d,a) \z\- d ~ a 


Xa,X ._ ^0 + a yA/a 


A + a^ 2 


a a A(d,a)\z\- d ~ a l {lzl<x} 


X a := I"- 1 


A + a a K a J 2 


a a A(d,a)\z\~ d ~ a l {lzl<1} 



Recall that p Q {x) := x d -(f) Q (x) for every Q e 3D and x e {y = (y, y d ) G B(Q, R) :y d > (pQ{y)}. 
We define for r\ , r<i > 

DQ(n,r 2 ) :={y e D :n> p Q (y) > 0, |y| < r 2 } . 

Lemma 3.4 There are constants 5q = 5o(R,M,A,a) G (0,ro), Cs = Cg(R, M, A, a) > and 
C 9 = C 9 (R,M,A,a) > suc/i that for every a G (0,M], A > 1, Q £ 3D and x G Dq(A _1 5o, A _1 r ) 
-urai/i x = 0, 

Gl?Q(2A^ 1 5o,A- 1 ro)') > C 8 A,5 D (x), (3.26) 

DQ(A-l5 ,A-lr ) / 



G D < C 9 A5 D (x) (3.27) 

D Q (A-lA ,A-lr ) 



and 



^(A-^cA-Vo) 



< C 9 A~%(cc). (3.1 



Proof. To derive the estimates in the lemma, it will be convenient to consider the scaled process 
XX^_2 t , which has the same distribution as X a ^ a 2> a ' X . The latter has infinitesimal generator 
A + a a X a - 2 A^ x 2 . 
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Without loss of generality, we assume Q = and let <fi : M. d 1 — ► R be the C 1 ' 1 -function satisfying 
0(0) = V(/»(0) = 0, ||V0[|oo < A, \Vc/){y) - V0(2)| < A|y - 2| and CSq be the corresponding 
coordinate system such that 

B(Q, R)nD = {(y, y d ) G B(0, 12) in CSq : y d > <f>{y)}. 

Note that, since D is a C ' open set with characteristics (12, A), for every A > 1, XD is a C ' 
open set with the same characteristics (12, A). Let 4>\(y) := </>(A _1 y) : R d_1 — > R. Then A satisfies 
0a(O) = V0 A (O) = 0, || V0a||oo < A, |V0 A (y)-V0 A (2O| < A|y - 2] and 

B(Q, R)n\D = {yG B(0, R) in CS Q : y d > 4>\{y)} for all A > 1 . 

We let p > be such that p/a and 1 < p < (2 A (3 — a)), and define 



pa(s/) 

h x (y) 

h \, P (y) 

D(X, ri ,r 2 ) 



PA(y)lB(0,4ro)nAD(y), 

h\(y) p = {p\(y)) P lB(0,4ro)nAD(y)) 

{y e AD : < p A (y) < r x and |y| < r 2 } . 



Since p\(y) < y/l + A 2 5xD(y) in view of (13, 2ft . we have < h\ < 12 < 1. It is easy to see that 
D(X,r\,r2) is contained in Dn 13(0, 12/4) for every r%,r2 < ro- Note that the (vector-valued) 
Lipschitz function V(j)\ is differentiable almost everywhere. So for a.e. y £ _B(0, 4ro) H AD, 

A/i A (y) = A(y d - <f> x (y)) = -A0 A (y) (3.29) 

and 

Ah x>p (y) = A(y d -<t> x (y))P 

= p(p- 1)(1 + |V0 A (y)| 2 )(p A (y)f- 2 -pipxwr-^m 
> p (p- l)(l + |V<My)| 2 )( PA (y)r 2 -p^^r^A^iu. 

Thus, since p € (1,2), we can choose a positive constant 5i = 5i(R, M, A,a) £ (0,ro), independent 
of A, so that there is c\ > such that 

&h x , p (y) > ci (p A (y)) p ^ 2 > for a.e. y e D(X,5i, r ) . (3.30) 

We divide the rest of the proof into three steps. 

Step 1: Constructing suitable superharmonic and subharmonic functions with respect to A + 
opiyoi— 2 /V?/ 2 _ Lgt, ip be a smooth positive function on R ' with bounded first and second order 
partial derivatives such that ip{y) = 2 p+1 |y| 2 /rQ for \y\ < rg/4 and 2 P+1 < ij)(y) < 2 P+2 for 
\y\ > ro/2. Now we consider 

ui,x(y) ■= h x {y) + hx, p (y) 

and 

«2,A(y) := h x (y) + tp(y) - h x , p {y). 
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Observe that since < h\ < 1 and p > 1, both ui A and U2,\ are non-negative. By Taylor's 
expansion with remainder of order 2, 



A + a a X^A a d ^) i>{y) < |A^(y)| + M 



Ad 7 A V(y) <c 2 (a,M) <oo. 



(3.31) 



Note that the constant c 2 above is independent of A. Moreover, since A > 1, p > a/2 and p ^ a, 
by (|3.16p and f|3. 1T|) there exist C3 = ^(R, A) > and 82 = 82(11, A) G (0, 5i] independent of A such 
that 

K!x hx,p(v) > -csA p " a for y G D(A, J 2 , r ). 

Thus by using (|3.30p . the fact that p < 2 and the inequality above, and by choosing 82 smaller if 
necessary, we get 



(A + a^A^-^A^ 2 ) h x>p (y) > c lPx (yf^ 2 - M a c^ p ~ a ^ a -^ 



(3.32) 



>c lP M p - 2 -M a cz> c ±p x (yy- 



for a.e. y G D(X,8 2 ,r ). Furthermore by (l3TT6l) - (l3TT8| ) and (I3T291) . there exist c 4 = c 4 (M) > and 
83 G (0, 82) independent of A > 1 such that for a.e. y G D(\, 83, tq) 



A + a^X^A^) h x (y) 

< C4 + A (l-«)++(a-2) pA(y) (l-«)A0 + 1{q=i} A -l| log(/0A(y)/A) |) (3.33) 

< C 4 (l + A (l-«) + +(«-2) pA(y )(l-a)A0 + e -l + 1{q=1} I log PA ( 2/ )|) . 



Thus by (|3,31|) - (|3.33|) and the fact that p < 2 A (3 — a), there exists 84 G (0, £3) independent of 
A > 1 such that 

(A + a-A^A^ 2 ) u 2 ,x(y) < c 2 + c 4 (2 + I log p x (y)\ + px^W ) - f p A (y) p " 2 < -1 (3.34) 

for a.e. y G D(X, 84, ro). 

On the other hand, we have from (|3.17p and (|3.18p . 



, + a a X^A a d ; x 2 )h x (y) > 

> 



■II A^aIIoc - c 5 M a (X^+^ + A" 1 log A + A" 1 ) logp A (y)|) 
■IIA^aIIoo - c 5 M a (l + e" 1 + I logp A (y)|) 



for a.e. y G D(X, 84, ro). Combining the inequality above with (|3.32p . by choosing 84 smaller if 
necessary, we have for a.e. y G D(A, #4,ro), 

(A + a^A^A^ 2 ) u hX (y) > -HA0IU - c 5 M a (2 + | logp A (y)|) + |p A (y) p - 2 > 0. (3.35) 

Step 2: Translating super- / sub-harmonic functions into super- / 'sub-martingale properties for X aA( ™ 2>//a > A . 
For notational convenience, we let 



X 



X 



lX (c-2)/ a X 



and r, 



a. A 

u 



lA («-2)/a X 
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We claim that the estimates f|3.34|) and (|3.35|) imply that 



t l— ^ u 2,\ ( ^ a '~ a ,A ) is a bounded supermartingale, (3.36) 



Mr ' 



E, 



~a,X 

T D(\,S 4 ,r ) 



< Px(x), (3.37) 



and 



t l— * u i,x( ^ a,X ~a,\ ) is a bounded submartingale. (3.38) 



tAr ' 



Observe that if v is a bounded C 2 -function on ]R rf with bounded second order partial derivatives, 
then by Ito's formula and the Levy system (|3,23|) , 



M? = v(X?> X )-v 



(A^' A ) - J (A + a Q A (a - 2) A^ A ) v{X^' X )ds (3.39) 



is a martingale (see the proof of Proposition 4.1 in [6] for the derivation of a similar assertion). If 
the functions U2,x and u\^\ were C 2 with bounded second order partial derivatives, then the claims 
(HEMP , ([337]) and ([3T38]) would just follow from ([3T39]) and the estimates (pT34"|) and (|3735]) . However 
they are not C 2 since D is C 1,1 and they are truncated on the outside of 5(0, 4ro) PI AD. So we 
will use a modifier. Let g be a non-negative smooth function with compact support in R rf whose 
value only depends on \x\ such that g{x) = for |x| > 1 and J Rd g(x)dx = 1. For > 1, define 
S fc (x) = 2 kd g(2 k x). Set 

u ! a( z ) := (9k * u i,x)( z ) '■= / 9k(y)ui,\{z - y)dy, i = 1,2. 

As 



A + a"A(- 2 ) A^ 2 ) «g = 9 , * (A + a^" 2 ) A« /2 ) u hX for i = 1, 2, 
we have by ff3T3"^ and (f!P5D that 

;A + a^- 2 )A^).g>0 and (a + a^" 2 ) A^ 2 ) < -1 



on D k (X,5 A ,r ) := {y : 6 4 - 2 k > Px (y) > 2 k and \y\ < r - 2 

Since i = 1,2, are bounded smooth functions on IR d with bounded first and second order 
partial derivatives, it follows from (|3.39|) that 

and 



CAT D fc (A,,5 4 ,r ) 



(A,<5 4 ,r ) * s a P os hive supermartingale 



i 1 — ► «i if^'-o.A ) i s a bounded submartingale. 

(k) 

Since for i = 1, 2, u^a is bounded and continuous, «W converges uniformly to u^a- Thus 

fc (A,<5 4 ,r ) is a Positive supermartingale (3.40) 



and 



i ► ^i,a ( X a,X ~ a .\ ) is a bounded submartingale. 

tAT D k (\,S 4 ,r ) 
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Since Dk(X, #4, ro) increases to D(X, <!>4,ro), we conclude that (|3.36p and fj3.38|) hold. Moreover, for 
each fixed k > 1 and t > 0, we have from (|3.40p that 



E x 



tAT D fe (A,« 4 ,r ) 



+ iAr 



a, A 

D fc (A,<5 4 ,r ) 



~a,A 

r D(X,6 i ,r ) 



< '"2,a( 3; )- Since 



Since U2 y x > 0, by first letting k — > 00 and then i — > 00, we get E x 
£ = 0, tp(x) = and so U2\(x) < P\(x). This proves (|3.37p . 

S'tep 3: Deriving the desired exit distribution estimates by utilizing the super- /sub-martingale prop- 
erty. Since ifi > 2 P+1 on \y\ > ro and ip(x) = 0, we have by (|3.36j) . 

P\(x) > u 2 ,x(x) 

^a,A 



> E, 



/ v a,X 
U 2 ,X X' x 



; x 



D(A,<5 4 ,r )/ 'D(\,S 4 ,r ) 
a.X 



\ T D(X,S 4 ,r ) 

We also have from (|3.38l) 



G (XD)\D(X,oo,r ) 
G (A J D)\£»(A,oo,r ) 



Px{x) < px{x) + px{x) p = u liX (x) < E 
Combining the two displays above, we get 

G D(X,oo,r ) 



v a,X 



'D(A,« 4 ,r ) 



< ( X a J a \ G AD 

T D(A,<5 4 ,r ) 



X ^-~a,A 



T D(A,<5 4 ,r ) 

(/.A 



(3.41) 



G AD 



T D(A,<5 4 ,7- ) 



^a,A 

X I ^~a,A 



T D(A,<5 4 ,r ) 



G (XD)\D(X,oo,r ) 



2P+ 1 - 3 
>— — - -Pa(x')- 



By 



2 (2P+i _ 1) 

GD(A,cx),ro)\D(A,25 4 ,r ) 



=E fl! 



<E, 



a, A 

~a.A 

T D(A,i 4 ,r ) 
~a , A 

T D(A,« 4 ,r ) 



(3.42) 







tA( Q " 2 )/») a ^(d,a) 



~a , A 

T D(A,« 4 ,r- ) 



D(A,oo,ro)\D(A,25 4 ,r ) |X s a,A - y\ d + a 

a a X a ~ 2 A{d, a 



L {\XP*-y\<X} 



dyds 



D(A,oo,ro)\D(A,25 4 ,r ) \Xs' X ~ y\ d+a 

y\- d - a dy) E x 



dyds 



<c 6 A(d,a)a a X a - 2 
<c 7 A(d,a)a a X a ~ 2 
<c 8 E x 



D(A,oo,ro)\D(A,2<5 4 ,r ) 



-a, A 

r D(A,<5 4 ,r ) 



D(X,2S 4 ,r )\D(X,35 4 /2,r ) 



\yl d ~ a dy E, 



~a,A 

r D(A,<5 4 ,r ) 



— a, A 

T D(X,S 4 ,r ) 



--c 8 ¥ x [X^ x x 

T D(X,S 4 ,r ) 



a a X a - 2 A(d,a) 

D(X,28 4 ,r )\D(X,35 4 /2,r ) \X^' X - y\ d +<* ^-yKV 

G D{X, 25 4 , r )\D(X, 3<5 4 /2, r )) . 



dyds 
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Thus from ([33T>([332l) 



G D(\,254,r ) ) > cgpx(x). 

T D(\,S 4 ,r ) 



X 



(3.43) 



Recall that < h x , p < 1. If \y\ > r /2, then ip(y) > 2 P+1 , we have 

u 2 ,x(y) = i>{y) + h x (y) - h x>p (y) > i>{y) - h Kp (y) > 2? > 1 for y G B(0, r /2) c . 

Furthermore, for y G B(0, 4ro) such that £4 < p\(y) < 4ro, 

«2,a(j/) = 4>(y) + h\(y) - hx >p (y) > p\(y) - p\(y) p > c w , 

where cio G (0, 1) depends on £4 and R. By using the last two observations, it holds that U2,a > 
cio > on (XD) \ D(X,54,ro). Therefore, by (|3.36p we get 



P\{x) > u 2 ,\(x) > E a 



\ra,X 



'D(A,« 4 ,r ) 



> cxoP, ( X a J a \ 



e\D). 



'-D(A,<5 4 ,r ) 



Since the process {X(X^_ 2t — X^),t > 0} under P x has the same distribution as {X 
Xq A( ^ ' X ,t > 0} under F\ x , we have from (|3.43[) that for x G -Dq(A _1 (54, A _1 ro) 



X?a G -D 

D Q (A-l« 4 ,A-lr ) 



(3.44) 



aA (a-2)/a x 



> 



, , €D Q (2X- 1 5 4 ,X- 1 r )\D Q (X- 1 5 4 ,X- 1 r )) 
G J D(A,2 ( 5 4 ,ro)\ J D(A,5 4 ,ro) 



^a,A 



"Ax I AT_ a A 

r D(A,44,r ) 

> c 9 pa(Ax) > c u 5 X d(Xx) = c n X5 D (x), 



and, from (|3.44j) 



AT"a G D 

D Q (A-li 4 ,A-lr ) 



I X^a G AD 

T D(X,S 4 ,r ) 



< c 12 p\(Xx) < ci 3 5 X d(Xx) = c 13 X5 D (x). 



Finally by (13351) and (13371) 

E 



r D Q (A- 1 5 4 ,A- 1 ro) 



G 



a.l 



r> Q (A- 1 5 4 ,A- 1 r ) 



£)q(A- 1 < 54 1 A-Vo) 



[x,y)dy 



A 



d-2 



^D(t5 4 2 , ) C)' A ( Ax ' Ay ) dy 
D Q (A-i5 4 ,A-ir ) V ° ; 



lD(\,S 4 ,r ) V ° ; 



A" 2 E Ax 



-a, A 



'D(A,54,r )_ 

< A~ 2 pa(Ax) < ci 4 A~ 2 5ad(Ax) = ci 4 A _1 (5z)(x). 
This completes the proof by taking So = 64, Cs = en, and C9 = max{ci3, C14}. 



□ 



We now derive exit distribution estimates for the process X a from those for X a in Lemma 137 
Recall that r = R/(4y/T+ A 2 ). 
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Lemma 3.5 There are constants 5q = 6q(R, M, A, a) G (0, ro), Cg = Cg(R, M, A, a) > and C\q = 
C 10 (R,M,A,a) > such that for every a G (0,M], A > 1, Q G dD and x G Dq(\~ 1 5 , A _1 r ) with 
x = 0, 

F x (x% eD Q (2X- 1 5o,\^ 1 r )] > C 8 \5 D (x), (3.45) 



and 



x , G £> J < Ci A5 D (x) 

V D Q (A-l« ,A-lr ) y 

< doA- 1 ^^). 



(3.46) 
(3.47) 



Proof. Without loss of generality, we assume Q = and let : R d_1 -» R be the C^-function 
satisfying 0(0) = V</>(0) = 0, ||V0||oo < A, | V<f)(y) — V(f)(z)\ < A\y— z\ and CSq be the corresponding 
coordinate system such that 

B(Q, R) n D = {(y, y d ) G CSq n 5(0, i?) : y d > cj>(y)}- 

Let (5o, Cg and Cg be the constants from the statement of Lemma l3.4i Since diam(Dg(A~ 1 5o, A _1 ro)) < 
i, we have that 



\x-y\ d a l{|x-y|<i} = \x-y\ d ~ 



for all x,y G Dq(\ 1 5 ,X V ). 



Let 



j(x) -a^a)^!-^)^,,^!}. 

Note that J^ d j(x)dx < oo. Thus we can write X" = Xf + where is a compound Poisson 
process with the Levy density j(x), independent of X^ . Since the jump size of Z a is greater than 
or equal to 1 and diam(I?Q(A _1 5oi A _1 ro)) < \, we see from ([3.28P that 



E z L r D Q (A-i5 ,A-ir () ) 

Moreover we have from (|3.26p that 



<E, 



r D Q (A-i<5o,A-ir ) 



^CgX-^Dix). 



Pxfc , , G D Q (2A- 1 (5o,A- 1 ro)) 
=W x (x2a 1 G D Q (2A- 1 5o,A- 1 r )) > C 8 \5 D (a 

We recall the notations from the proof of the previous lemma: 

P\{x) ■= yd ~ ^A" 1 ?/), 
D(X,r 1 ,r 2 ) := {y G CSq : n > p x (y) > 0, \y\ < r 2 }, 



x a ' A = x aX(a ~ 2)/0 '' x 



and r, 



a,A ^aA< a - 2 )/ Q ,A 
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Let a (A) := aA (a_2)/a , which is no larger than M. By ([L~4"|) . 



E, 



a(A) 
a(A) 

D(A,S ,r ) 
a(A) 

C(A,5 ,r ) 



G (A J D)\ J D(A,2 < 5 ,2r ) 

a(A) a ^(d,a) 



J (XD)\D(X,2S ,2r ) \Xs W - y\ d + a 



dyds 



< Cl A(d,a)(a(X)) a 

< c 2 A(d,a)(a(X)) a 

< c 3 E x 



\y\- d ~ a dy E x 



(\D)\D(\,2S ,2r ) 



a(A) 
r D(\,8 ,r ) 



\yr d ~ a dy E x 



D(\,28o,r )\D(\,36o/2,r ) 



a(A) 
r D(\5 ,r ) 



a(A) 

D(A,«o^o) 



a(A) a ^4(d,a) 

D(A,25o,T )\Z)(A,3<5o/2,ro) |Xs (A) - y\ d+a 



dyds 



= c 3 ¥ x (Vg>> G L»(A, 25 , ro) \ D(X, 35 /2, r )) 
Thus by the above inequality and (|3.44p . we have 

r a(X) 



G AD 



O(A,5 ,r ) 
,a(A) 



a(A) 

r r>(A,« ,r ) 



< c 3 P. ( X a lfi 



XTd) e (A-D) \ £>(A, 2«5 , 2r ) ) + P. ( X!; A A G D(A, 2<5 , 2r ) 



'D(A,« ,r ) 



O(A,«5 ,r ) 



G D(X, 25 Q , r ) \ D(A, 3S /2, r ) ) + P. ( G AD 



c 3 Px ( *!; A 



T D(A,<5 ,r ) 



G D(X, 2S , r ) \ D(X, 35 /2, r ) ) + ( X% x € XD 

T D(\,S ,r ) 



D(X,5 ,r ) 
Ja,A 



< (c 3 + 1)F X ( X!; A A G XD J < c 4 p A (x). 

\ T D(\,S ,r ) / 

Since (AX? 2t ,i > 0) is the independent sum of a Brownian motion and a symmetric a-stable 
process on R rf with weight a(A), we have from (|3.49j) that for x G -Dq(A _1a o, A _1 ro) 



(3.48) 



(3.49) 



\[x% eD 

DgCA-iSo.A-lro) 



„a(A) 

Ax I ^ a(A) 



G XD 



'D(\,S ,r ) 

< Cipx(Xx) < c 5 (5ad(Ax) = c 5 X5 D (x). 



The proof is finished by taking Cio = max{Cg, C5}. 



□ 



4 Boundary Harnack principle 

In this section, we give the proof of the boundary Harnack principle for the independent sum of 
a Brownian motion and a symmetric stable process. We first prove the Carleson estimate for the 
independent sum of a Brownian motion and a symmetric stable process on Lipschitz open sets. 



25 



We recall that an open set D in M. d (when d > 2) is said to be a Lipschitz open set if there 
exist a localization radius R\ > and a constant Ai > such that for every Q G dD, there exist 
a Lipschitz function 4> = 4>q : W 1 ^ 1 — > R satisfying </>(0) = 0, |</>(x) — </>(y)| < Ai|x — y\, and an 
orthonormal coordinate system CSq: y = (y±, . . . , yd-i, yd) ='■ (y, Vd) with its origin at Q such 
that 

B(Q, R l )DD = {y = (y, y d ) G B(0, R x ) in CS Q : y d > <f>(y)}. 

The pair (i?i,Ai) is called the characteristics of the Lipschitz open set D. Note that a Lipschitz 
open set can be unbounded and disconnected. For Lipschitz open set D and every Q G dD and 
x G B(Q, Rx)r\D, we define 

p Q (x) := x d - 4> Q {x) , 

where (x,x d ) is the coordinates of x in CSq. 

We recall that X" = X® + aYt is a Levy process with characteristic exponent <£ a (x) = \x\ 2 + 
a a |x| a . This process may be obtained by subordinating a d-dimensional Brownian motion W = 
(Wt, t > 0) by an independent subordinator T t a := t + a 2 T t where T = (T t , t > 0) is an a/2- 
stable subordinator. More precisely, the processes X" and Wt? have the same distribution. Note 
that the Laplace exponent corresponding to T a is equal to 4> a {\) = A + a a \ a / 2 . Let M a /2(t) := 
Sn^=o( — l)"i na//2 /r(l + na/2). It follows by a straightforward integration that 

/■OO -I 

which shows that the potential density u a of the subordinator T a is given by 

u a (t)=M l _ a/2 (a 2a ^ 2 ~ a h). 

Since, for any a > 0, <p a is a complete Bernstein function, we know that u a (-) is a completely 
monotone function. In particular, u a {-) is a decreasing function. Since u a (t) = u 1 (a 2a// ( 2_a )t), we 
know that a \— > u a (t) is a decreasing function. Therefore, if < ai < ci2> then u ai {t) > u a2 (t) for 
all t > 0. We will need this fact in the proof of next lemma. 

Lemma 4.1 Let D C M d be a Lipschitz open set with the characteristics (Ri, Ai). There exists a 
constant 5 = S(Ri, Ai, M) > suc/i i/taf /or a// a G [0, M] and Q G x £ D with pq(x) < R\/2, 

F x (X« {x) eD c )>6, 

where r(x) := r a Dr]B{x ^ pQ(x)) = M{t > : X? $ D n B(x, 2p g (s))}. 

Proof. Clearly, 

F x (X a r{x) G D c ) > F x (X a r[x) G D c n 2p (x))) > P, (x« (x) G (91? n 2p Q (x))) . 

Let D x := DPiB(x, 2pq{x)) and be the subprocess of Brownian motion W killed upon leaving 
D x . The process Z a defined by Zf := W Dx (Tf), where T t a is an independent subordinator described 
in the paragraph before the statement of the lemma, is called a subordinate killed Brownian motion 
in D x . We will use £ to denote the lifetime of Z a . It is known from [39] that 

F x (x? (x) G dD H B(x, 2p Q (x))) > P. G 9-D n 2p g (x))) 

= E x [ U a (T Dx ); W fDx edDH B(x, 2p Q {x))} . 
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Here and below, tjj := inf{t > : Wt ^ U} is the exit time of W from U. Denote C x := 
dDPiB(x,2p Q (x)). Then 

E x [u a (? Dx );Wr Dx GC x ] > E x [u a (r Dx );Wr Dx GC x ,r Dx <t] (4.1) 

> u a {t)f x [Wr Dx GC x ,r Dx <t] 

> u a (t) (¥ x (W 7dx G C x ) - F x (t Dx > t)) , 

where t > will be chosen later. 

Since D is a Lipschitz open set with characteristics (R±,Ai), there exist r\ = rj(A,Ri) > and 
a cone 

C:={y = ( yi ,...,y d )£R d :y d < 0, (y? + • • • + 2/ii) 172 < v\Vd\} (4.2) 

such that for every z G <9-D, there is a cone with vertex z, isometric to C, satisfying C z CiB(Q, Ri) C 
D c . Then by the scaling property of W and symmetry considerations, we have 

F x (w 7Dx eC x ) > P x {W 7B(Xt2pQ{x)) e dB(x,2 PQ (x))nC { ^ Q(B)) ) 
> p o(^ s(0i2) g di?(o,2) n(C + (o,-i))), 

which is strictly positive. Hence we can conclude that there exists c\ = c\{D) > such that 

f x (W 7Dx G C x ) > a . (4.3) 

Next, 

F x (r Dx >t)< ^1 < gi^gg^ < C2 ^M! < C2 M , (4.4) 
for some constant C2 > 0. By using (14. 3p and ()4.4j) in (j4.1j) . we obtain that 

Ex [u a (T Dx )- Wt Dx G C x ] > U a (t) (cy - C 2 ^pj . 

Now choose t = Ai) > large enough so that ci - c 2j Rf/t > ci/2. Then 

E x [«(tdJ; Wr Dx G Cj > ciu a (i)/2 > Cl u M (t)/2 =: 5. 

The lemma is thus proved. □ 

Suppose that D is an open set and that U and V are bounded open sets with V C V C U and 
If u vanishes continuously on D c n £/, then by a finite covering argument, it is easy to 
see that u is bounded in an open neighborhood of dD n V. 

Lemma 4.2 Let D be an open set and U and V be bounded open sets with V C V C U and 
D fl V 7^ 0. Suppose u is a nonnegative function in M. d that is harmonic in D C\U with respect to 
X a and vanishes continuously on D C (~)U. Then u is regular harmonic in D(~)V with respect to X a , 
i.e., 

u(x) = E x \u(X* a )] for all x £ DDV . (4.5) 
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Proof. For n > 1, let B n = {y E D n V : 5d{u) > V n }- Then for large n, B n is an non-empty 
open subset of D n V whose closure is contained m. D P\U . Since u is harmonic in D n U with 
respect to X a , for x £ D C)V and n large enough so that x E B n , we have that 



u(x) = E x u(X% ) 



E 3 



u 



, A r" ) > r B n — r £>nv 



Hence 



u(x)-E x u(X?« nv ) 



(4.6) 



< E, 



+ E, 



Since lim n ^oo = T^ nV P x -a.s., the second term in (|4.6p converges to E x [ti(X"a ); A] where 
A := n^ =1 {r| n < r£ nv }. Note that 

E 3D n V on A. 

Hence u =0 on A, as u is assumed to vanish on D c n £/. Consequently 

^ 'nnl/' 



lim E T 



0. 



For the first term in (|4,6p . note that 5d{X^ b ) < 1/n on { r |? n < T hn\A- Therefore, by the 
assumption that u vanishes continuously on D C D U, one has lim n ^ oc u(X% g ) = 0. Moreover, since 
u vanishes continuously on (dD) n U, there is no > 1 so that u is bounded in D n V \ S no . So by 
the bounded convergence theorem we have 



lim E, 



u 



( X r B J; T B„ < T iW 



0. 



This proves the lemma. 



□ 



Proof of Proposition ll.3i We know from the parabolic Harnack inequality from Theorem 6.7 
of [20] that Harnack inequality holds for the process X := X 1 . That is, there exists a constant 
ci = c\(a, M) > such that for any r E (0, M a ^ 2 ~ a ^], xq E M d and any function v > harmonic 
in B(xq, r) with respect to X, we have 



■u(x) < c\v(y) for all x, y E B(xq, -). 



(4.7) 



Now the proposition is an easy consequence of (|4.7j) . In fact, note that for any a E (0, M], X a has 
the same distribution as XX \-2 t , where A = a a ^ a ~ 2 ^ > M a ^ a ~ 2 \ Consequently, if u is harmonic 
in B(xo,r) with respect to X a where r E (0,1], then v(x) := u(Xx) is harmonic in B(X~ 1 xo, A r) 
with respect to X and A~V < M a /( 2 " a ). So by (|477]) 



That is, 



u(Ax) = f(x) < ciu(y) = ci?/(A?/) for all x,y £ B(X xq,X 1 r/2). 



u(x) < c\u{y) for all x,y £ B(xo,r /2). 



□ 
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Theorem 4.3 (Carleson estimate) Let D C R d be a Lipschitz open set with the characteristics 
(i?i,Ai). Then there exists a positive constant A = A(a,A\,Ri,M) such that for a G (0, M], 
Q G dD, < r < R\/2, and any nonnegative function u in W 1 that is harmonic in D n B(Q,r) 
with respect to X a and vanishes continuously on D c n B(Q,r), we have 

u(x)<Au(x ) for x e Df]B(Q,r/2). (4.8) 

where xq G D n -E>(<2, r) wii/i pq(xo) = r/2. 

Proof. Fix a £ (0, M]. Since I? is Lipschitz and r < R±/2, by the uniform Harnack principle in 
Proposition 11.31 and a standard chain argument, it suffices to prove (|4.8p for i £ flfl B(Q, r /12) 
and xq = Q. Without loss of generality, we may assume that u(xq) = 1. In this proof, constants 
5, P, rj and Cj's are always independent of r and a. 
Choose 0<7<a/(d + a) and let 

B = D n B(x, 2p Q (x)) , 5i = £?(x, r 1 '^ p Q { X y) . 

Further, set 

B 2 = B(x , pq(x )/3) , B 3 = B(xq, 2p Q (x )/3) 

and 

r = mi{t > : X? i B } , r 2 = inf{i > : X? $ B 2 }. 
By Lemma l4.1| there exists 5 = 6(R±,Ai, M) > such that 

F X (X? £D C )>5, x G B(Q,r/A) . (4.9) 

By the uniform Harnack principle in Proposition 11.31 and a chain argument, there exists (3 such that 

u(x) < (p Q (x)/r)- l3 u(xo) , x G DD B(Q,r/A) . (4.10) 

In view of Lemma 14.2} u is regular harmonic in Bq with respect to X a . So 

u{x) =E x [u(X? );X? G B 1 ) + E x ); X% $ B x \ for x G B(Q,r/A). (4.11) 

We first show that there exists r] > such that 

E4«(jq );X^ ^Si] <u(x ) if xe DnB(Q, r/12) with p Q (x)<r l r. (4.12) 

Let r/o := 2~ 2( - d+a ^ d , then for />q(x) < r/ r, 

(p Q (x)) d /^ < 1/4 and 2 PQ (x) < r 1 ^ p Q ( X y - 2p Q (x). 

Thus if x G -D n B(Q,r/12) with pq(x) < r/or, then |x — y| < 2|z — y| for z G -Bq; V & B\. Thus we 
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have by (jl,4p and Lemma 13. II 

E x [u(X? );X? iB x ] (4.13) 



= A(d,a) / G% (x,z) / a a \z-y\- d - a u{y)dydz 

JB J\y-x\>r 1 -''p Q (x)~< 

< 2 d+a A(d,a) [ G a B() (x,z)dz [ a a \x-y\- d - a u(y)dy 

JB J\y-x\>r 1 -~fp Q (x)"< 

< d+a a r^_, _ , ,,l / „<*\„ _ n ,\-d-a 



2 d+a A(d, a)E x [T B ^ 2pQ{x)) ] / a a \x - y\- d - a u(y) dy 

\y-x\>r 1 -r p Q (xp 



< 2 d+a A(d,a)c lPQ (x) 2 ( [ 

\J\y~x\yr 1 - 



a a \x - y\ d a u(y) dy 

">p Q {x)-<,\y-x \>2p Q {x )/3 



+ / a a \x-y\- d - a u(y)dy 



' B2 J\y-x \>2p Q (x )/3 

a a \x -y\- d - a u(y)dy 



= ■ c 2PQ {xY{h+h). 

On the other hand, for z £ B 2 and y ^ B3, we have \z — y\ < \z — xo\ + \xo — y\ < pq(xo)/3+\xq — y\ < 
2 |^o — y\- we have again by (|1-4|) and Lemma 13. II 

u(x ) > E X0 [ U (X^),X^ 2 ^B 3 ] (4.14) 

> A(d,a) f G a B2 (x ,z) I a a \z-y\- d - a u(y)dydz 

JB 2 J\y-x \>2p Q (xn)/3 

> 2~ d - a A{d,a) [ G a B2 (x ,z)dz [ 

JB 2 J\y-xo\>2p Q (xo)/3 

> 2- d - a A(d,a)c 3 ( PQ (x )/3) 2 [ a a \x -y\- d - a u(y)dy 
= c 4 pq(x ) 2 / a a \x -y\~ d ~ a u(y)dy . 

J\y-x \>2p Q (x )/3 

Suppose now that \y — x\ > r 1_7 pQ(x) 7 and x G B(Q,r/A). Then 

\y - xo\ < \y - x\ + r < \y - x\ + r 1 p Q (x)~ 1 \y - x\ < 2r 7 p Q {x)~ 1 \y - x\. 

Therefore 

h= f a a \x - y\- d - a u(y) dy (4.15) 

J^-x^r 1 —' p Q (x)i ,\y-x \>2p Q (x )/3 

< [ (2- 1 (p Q (x)/rr)- d - a a a \y-x \- d - a u(y)dy 

J\y-x \>2p Q (x )/3 



= 2 d+a (p Q {x)/r)-^ d+a n a a \y-x \- d - a u(y)dy 

Av- x o\>2pq(x )/3 

< 2 d+a (p Q (x)/r)-^ d+ ^C4 1 p Q (x )- 2 u(xo) 
= c 5 (p Q (x)/r)-^ d +^p Q (x )- 2 u(xo), 



where the last inequality is due to (|4.14p . 
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If \y — xq\ < 2/3q(xo)/3, then \y — x\ > \xo — Q\ — \x — Q\ — \y — xq\ > pq(xq)/Q. This together 
with the uniform Harnack principle in Proposition 11.31 implies that 

h = I a a \x-y\- d - a u(y)dy (4.16) 

J\y~ x o\<2p Q (x )/3 

< c 6 / a a \x - y\~ d - a u(x ) dy 

^l«/-a;ol<2p(3(xo)/3 

< c 6 u(x ) a a \x - y\~ d ~ a dy = c 7 a a p Q (x )~ a u(x ) . 

J\y-x\>p Q (x )/6 

Combining (|4.13p - (|4.16p we obtain 

E.Ki;);!;^!] (4.17) 

< c 2 p Q (x) 2 (c b {p Q {x)/r)-^ d ^ p Q (x r 2 u(xo) + c 7 a a p Q (x )- a u(x i 

< c 8 u(x ) (p Q {xf{p Q {x)/r)-^p Q {x,y 2 + a a p Q {xfp Q {x») 

< c 9 u(x )((p Q (x)/r) 2 ^ d+ ^ +a a p Q {xfr- a ) , 
where in the last inequality we used the fact that pq{xq) = r/2. Choose now r/ G (0,170) so that 

c 9 (y-T( d + Q ) + r} 2 M a \ < 1 . 
Then for x £ D D B(Q, r /12) with pq(x) < rjr, we have by (|4.17|) 

E x X? i B 1 ] < c 9 u(x ) (if-^^ + r, 2 r 2 - a M^ 

< C 9 (rf~^ d+a ) + V 2 M a ^ U(X ) < U(X ) • 

We now prove the Carleson estimate (|4.8p for x G Df]B(Q, r/12) by a method of contradiction. 
Recall that u(xq) = 1. Suppose that there exists x\ G D n B(x,r/12) such that u(x\) > K > 
r]~P V (1 + 5~ 1 ), where if is a constant to be specified later. By (|4.10p and the assumption 
u(x\) > K > r]~P, we have (pQ(xi)/r)~P > u(x\) > K > r}~@, and hence Pq{x\) < rjr. Let Bq, 
B\ and tq be now defined with respect to the point x\ instead of x. Then by (|4.1ip . (|4.12p and 
K > l + 5~\ 

K < u( Xl ) < E X1 [u{X%);X^ GBjJ+1, 

and hence 

E X1 [u{X a To )-X a T0 G B,} > u(xi) - 1 >—«(!!) . 

In the last inequality of the display above we used the assumption that u(x\) > K > 1 + d^ 1 . If 
K > 2 f3 /~ / , then D c f)Bi C D c n B(Q, r). By using the assumption that u = on D c n 5(Q, r), we 
get from (|4.9p 



E X1 [«(X° G Si] = K xl [u(K)>K £B 1 DD]< V X (X? G D) supu < (1 - 8) supu. 

Bi B 1 

Therefore, sup Bl u > u(xi)/((l + <5)(1 — 5)), i.e., there exists a point 22 G D such that 



ki-a^l^r 1 7 / oq(xi) 7 and u(x 2 ) > u(xi) > — 

1 — o z 1 



if. 
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By induction, if x k G D n B(Q,r /12) with ii(xfc) > 2f/(l - 5 2 ) fc_1 for > 2, then there exists 
x k+ i G 25 such that 

\xk ~ %k+i\ < r 1 ~ 7 pQ(x k )' y and u(x k+l ) > - ^ ^ u{x k ) > ZTp)k K ■ ( 4-18 ) 
From (|4~T0j) and (flTHj) it follows that PQ{x k )/r < (1 - <5 2 ) (fe ~ 1)//3 K- 1//3 , for every k > 1. Therefore, 

fc— 1 oo 

3=1 w 3=1 

oo 

< — + r 1 " 7 V(l - ( 5 2 )(j- 1 )T/^ET" 7//3 r 7 
12 3=1 

oo 

_ — + r l ~~t r ~i / P ^(1 - S 2 )^"'/ 13 

12 3=0 

• +ri ^-7//3. 



12 1- (l-J2)7//3 - 

Chooser = 7 ? V(l + J- 1 )Vl2^ 7 (l-(l- ( 5 2 ) 7 /^)-^/ 7 . Then K~ 7 ^ (1 — (1 — d 2 ) 7 ^) -1 < 1/12, and 
hence x& G DC\B(Q, r/6) for every > 1. Since lim^oo u(x k ) = +oo, this contradicts the fact that 
u is bounded on B(Q,r/2). This contradiction shows that u(x) < K for every x G Df]B(Q,r/12). 
This completes the proof of the theorem. □ 

Proof of Theorem ll.4l . Since D is a C 1 ' 1 open set and r < R, by the uniform Harnack principle in 
Proposition 1 1.31 and a standard chain argument, it suffices to prove (|1.6p for x, y G Df]B(Q, rro/8). 
In this proof, constants r\ and Cj's are always independent of r and a. 

We recall that ro = and <5o G (0, ro) is the constant in the statement of Lemma 13.51 

For any r G (0, 22] and x G DnB(Q, rro/8), let be the point Q x G 525 so that \x— Q x \ = <5d(x) 
and let xo := Q x + |(^ — Qx)/\x — Q x \. We choose a C 1 ' 1 -function : M^ 1 — > K satisfying 
0(0) = V0(O) = 0, ||V0||oo < A, |V0(y) — V0(z)| < A|y — z|, and an orthonormal coordinate 
system CS with its origin at Q x such that 

B(Q X , R)nD = {y= (y, y d ) G 5(0, 2?) in CS:y d > </>($)}. 

In the coordinate system CS we have x = and xo = (0, r/8). For any &i, 62 > 0, we define 

D(b 1: b 2 ) := {y = (y, y d ) in CS : < y d - </>(y) < hr5 /8, \y\ < b 2 rr /8} . 

It is easy to see that 25(2,2) C D n B(Q,r/2). In fact, since ro < ^ and <5o < for every 
z G 25(2,2) 

\z-Q\<\Q-x\ + \x-Q x \ + \Q x - z\< r - + r - + \z d - <f>(z)\ + |0(S)| < J(l + <5 ) < ^. 

Thus if u is a nonnegative function on R rf that is harmonic in D n 2?(Q, r) with respect to X a and 
vanishes continuously in D c f]B(Q,r), then, by Lemma [4.21 u is regular harmonic in D(~]B(Q,r/2) 
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with respect to X a , hence also in D(2,2). Thus by the uniform Harnack principle in Proposition 
11.31 we have 



u(x) = E x 



u(X a T a ' 



>E X 



u(X a T a VXta GD(2,1) 



(4.19) 



> ci«(a;o)P* i} G D(2, 1)) > c 2 u{x )5 D {x)/r. 



In the last inequality above we have used (|3.45p . 

Let w = (0, rro/16). Then it is easy to see that there exists a constant rj = r](A,ro,5o) £ (0, 1) 
such that B(w,rjrr /16) € D(l, 1). By (fL4"j) and LemmaEEl 



u(w) > E w u(X" a );X? a <£D(2,2) 

y ' ~ [ y T D(1,1) J D(l,l) r V ' ' 

= A(d,a)a a [ G a D[11) (w,z) [ 

JD(IA) JR d 



u(y) 



> 



c 3 a a E w [r% 



B(w,r]rro/(16))_ 



p(2,2) \z-y 

U{y) dy 



d+a 



dydz 



,2) \w-y 



d+a 



> c A a a r 2 



u(y) 



,2) \w-y 



d+a 



dy. 



Hence by (^47]) . 



E x 



u[X% )]X% <£D(2,2) 

V -D(l,l) / 



A(d,a)a a I G a D{11) (x,z) 

D(l,l) ' JM d \D(2,2) \ z ~ V 

u{y) 



u(y) 



d+a 



dydz 



< c 5 a a E x [r a D{ )] / 



2) f -y 
u(y) 



d+a 



dy 



2)\w-y 



c 6 5 D (x) 
dy < u (w . 



On the other hand, by the uniform Harnack principle (Proposition ll.3j) and the Carleson esti- 
mate (Theorem 14. 3p . we have 



E, 



u[X% ) ; G £>(2,2)J < c 7 n(x )P, G #(2,2)) < c 8 u(x )fe(x)/r. 



In the last inequality above we have used (|3.46|) . Combining the two inequalities above, we get 



u{x) 



E, 



+E 2 



u[X% )]X% eD(2,2) 

\ 15(1,1) / 

nfc );X"a $ D(2,2) 

\ T D(1,1)/ £>{1,1) 



(4.20) 



< — Od\x)u(xq) H 

r C4 r 

< — (5d(x)(u(xo) + u(w)) 
r 

cio 



< 



5d(x)u(x ). 



In the last inequality above we have used the uniform Harnack principle (Proposition 1 1 ,3| ) . 
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Prom (|4.19p - (|4.20p . we have that for every x, y G D n B(Q, rr /8), 

< Cio <fo(x) 



«(y) c 2 <5r>(y)' 

which proves the theorem. □ 
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